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ABSTRACT 

The prospects for accomplishing X-ray polarization measurements appear to have grown in recent years after a more 
than 3 5 -year hiatus. Unfortunately, this long hiatus has brought with it some confusion over the statistical uncertainties 
associated with polarization measurements of astronomical sources. The heart of this confusion stems from a 
misunderstanding (or potential misunderstanding) of a standard figure of merit — ^the minimum detectable polarization 
(MDP) — that one of us introduced many years ago. We review the relevant statistics, and quantify the differences 
between the MDP and the uncertainty of an actual polarization measurement. We discuss the impUcations for future 
missions. 
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1. INTRODUCTION 

Only a few experiments have conducted X-ray polarimetry of cosmic sources. In several rocket observations, we 
measured X-ray polarization from the Crab Nebula\ Using the X-ray polarimeter on the Orbiting Solar Observatory 
(OSO-8), we confirmed this result with a 19-sigma detection (P = 19.2% ± 1.0%)^, thus proving the synchrotron origin 
of the X-ray emission from this object. Data from OSO-8 also provided upper limits to polarization of a handful of bright 
targets. The key challenge for broad-band X-ray polarimetry has always been to obtain an instrument that is highly 
responsive to polarized radiation yet has low net background. In the past decade the development of a class of 
polarimeters based on the photoelectric effect^ '* and placed at the focus of an X-ray telescope has produced a renewed 
interest. 

2. STATISTICS 

All past and currently planned X-ray polarimeters for asfronomical applications share the feature that ultimately the data 
may be described as biimed as a function of the azimuth (cp) around the direction of the line of sight to the source. 
Polarization manifests itself as modulation at (2(p). Often polarimeters are rotated around the line of sight to address 
potential systematic effects; thus one commonly uses the language that polarization will manifest itself as a modulation 
of the detected signal at twice the rotation frequency. Here both amplitude and phase information are important as the 
former is related to the degree of polarization and the latter to the position angle. 

In the presence of noise, detecting this type of signal is most often discussed in classical post-war textbooks in the 
context of modulated radio waves in the presence of noise^. In the special case where the noise is described by Poisson 
distribution and is the total number of counts, the probability p(a, (p) of measuring a particular modulation amplitude a 
and a phase p, given that the true amplitude and phase are ao and po, is as follows: 



, , Na 
p{a,(p) = ——e^ 
An 



(1) 



' martin.c.weisskopf@nasa.gov; phone +1(256)961-7798; 
NASA/MSFCA^62, 320 Sparkman Drive, Huntsville, AL 35805-2718 USA 



It should be readily apparent, since the modulation is positive definite, that there is always a probability of measuring an 
amplitude of modulation, even if the underlying source is unpolarized (ao=0). Equation (1) may be integrated 
analytically and inverted to determine the amplitude which has a given probability of being exceeded by chance for an 
impolarized (ao=0) source. For example, the amplitude, that has only a 1% probability of being exceeded by chance 
is 
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One is not directly interested in the amplitude of modulation but the degree of polarization that corresponds to. There are 
two steps in relating these. The first step is recognizing that there may be a contribution from a steady background to the 
detected signal and therefore to express the modulation as a fraction of the detected signal only: 
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Here Rs and Rb are the signal and background counting rates. The next step introduces the modulation factor, M, the 
amphtude of modulation for a 100% polarized input in the absence of background. The modulation factor may be a 
function of energy. In the typical photoelectric polarimeter, M will be small at the low energies where tracks are short 
and the angular distribution appears almost spherical, but will rises very quickly above these energies to values possibly 
approaching 0.5. We define the minimum detectable polarization (MDP) as (a/A/): 
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By this definition, the "minimum detectable polarization" is the degree of polarization corresponding to the amplitude of 
modulation that has only a 1% probability of being detected by chance. The MDP, however, is not the uncertainty of a 
polarization measurement. This uncertainty, at a certain level of confidence, is found by integrating Equation (1): 



C = ^^a da d(p p(a, (p) . 
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Here C is the desired confidence level (0 <C < 1) and the integral is performed over an appropriate sub-region of the 
parameters space (a, cp). To invert Equation (5), we begin by defining the Stokes amplitudes: 



ao.x = do cos (po, 
ao,y = ag sin 
ax= a cos (p = ao,x cos tj/, 
ay= a sin ^ = ao,y +Aa sin i// 

Here \|/ is a parametric variable that we shall vary between and 2;r and 
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The amplitude a and phase angle (p ona. given confidence contour may then be expressed terms of ao, (po, y/, and Aa 
=Aac as 

a-{al+ Aa^ + la^ Aa^ cos(^^ - q)^'^ . 
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Using the relationship a da d<p = dAa^ dAa y= Aa dAa d i//, Equation 5 may be expressed as 
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The integrals in Equation 9 may be evaluated and solved for Auc- 



Aflc=J--ln(l-C) 
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Equations 8 and 10 then may be used to draw the appropriate confidence contours. 

3. HOW MANY COUNTS ARE NECESSARY? 

We are now in a position to examine the answer to two questions: (1) How many counts are required to achieve a 
particular minimum detectable polarization? (2) How many counts are required to measure this same level of 
polarization to a given accuracy? To simphfy the discussion we set the background level to zero. (Allowing for the 
background we leave as an exercise for the reader.) With no background, Equation 4 simplifies and may be solved for 
the number of counts, Ni necessary to achieve a particular MDP: 
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For example, Ni = 1,565,369 for a modulation factor of 1/3 and an MDP of 0.01 (1%). Next we want to answer the 
question how many counts are required to measure the amplitude of modulation corresponding to this MDP to say 3, 4, 
5, P-sigma accuracy? We have determined the answer to this question in two ways. In the first, we made use of the 
system of equations (8) and (10) to determine the extremes in the amplitude of modulation on the 67% confidence 
contours. We then foimd the number of counts, N2, such that 
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Where P =1, 2, 3, 4, ... is an approximate proxy for the number of sigma of the measurement. We fmd that the ratio, 
N2/N1 scales as /f and with a coefficient that for all practical purposes is 'A as shown in Figure 1. As a cross-check we 
performed the following derivation: 

W^i) = 4.29/Va^; 
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Fmally, let a„„,c - «mm,c = JJ Aa(c,N2); then 
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which reproduces Figure 1 for ri=4.29/2 and C=0.67. 
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Figure 1 . The increase in required observing time {N2/N1) versus tie number of sigma (JJ) in order accurately to 
measure the polarization. 
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Figure 2. The 1-a uncertainty in the position angle versus the number of sigma (yff) in figure 1. 



For completeness, we show in Figure 2 an estimate of the uncertainty in the position angle based on '/2 the extremes on 
the 1-sigma contours. This estimator scales as 1/(3. As with the degree of polarization, the conclusion that a precise 
measurement of the position angle requires far more counts than those necessary to determine the MDP is inescapable. 

4, SUMMARY AND CONCLUSION 

Here we have reviewed the statistical description appropriate for an X-ray polarimeter that attempts to measure 
polarization by searching for a variation of the detected signal at twice the rotation frequency. Our analysis applies to the 
case where the signal is dominated by Poisson statistics and our examples ignore the background. We have explicitly 
shown that the commonly used figure of merit for such a polarimeter — i.e., minimum detectable polarization (MDP) — is 
not the uncertainty in the polarization measurement. In particular, we compared the MDP at 99% confidence with a 
nominal measurement at the P-sigma level: We have shown that the number of counts required to measure the 
polarization is larger than that required to establish the 99% MDP, by a factor of p'/4\ Thus, a proper measurement, even 
at the 3-sigma level, requires an integration time at least 2.25 times longer than that simply to establish that one might 
have detected polarization at the level set by the MDP9(,%. A more reasonable criterion may be that the measurement 
should be accurate to 4- or 5-sigma if one wants to discriminate amongst different theoretical models for the source of 
the polarization: This requires an increase in observing time by a factor of 4 or 6.25, respectively. The implication is that 
observing programs where the integration time is set by the MDP may significantly underestimate the time required for 
measurement of polarization. This is especially true in those cases where there are grounds (prior measurements, 
physics, and observations at other wavelengths) to expect a certain degree of polarization to be present. For example, the 
degree of accuracy of polarization measurements as a function of pulse phase for the Crab pulsar may be estimated from 
optical measurements. 
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